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We develop a method to compute shakeup effects on photoluminescence in a Wigner crystal from localized 
holes. Our method treats the lattice electrons and the tunneling electron on an equal footing, and uses a quantum- 
mechanical calculation of the collective modes that is realistic throughout the Brillouin zone. We find that shakeup 
produces a series of sidebands that may be identified with maxima in the collective mode density of states, and 
definitively distinguishes the crystal state from a liquid state. We also find a shift in the main luminescence peak, 
that is associated with lattice relaxation in the vicinity of a vacancy. 
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Sixty years ago, Wigner [jjj pointed out that an elec- 
tron gas will undergo a zero-temperature, quantum phase 
transition into a crystalline phase as the density is low- 
ered. Forty-five years later, the first convincing evidence 
of an electron crystal was presented for a system of elec- 
trons on a He surface [||. The electron densities attain- 
able in this fashion arc extremely low, however, making 
this an unattractive system for observing the quantum 
phase transition. Semiconductors are much more attrac- 
tive systems in this sense, because one has great control 
over the electron densities, through dopant concentra- 
tions. A particularly good candidate for observing the 
Wigner crystal (WC) is the two-dimensional electron gas 
(2DEG), as realized in modulation doped semiconduc- 
tors. Samples of this type are now available with such 
high quality that the electron groundstate is not neces- 
sarily dominated by disorder. The possibility of observ- 
ing the WC is further enhanced by the application of a 
strong perpendicular magnetic field, which quenches the 
kinetic energy, and allows the formation of a crystal state 
at higher densities (for which disorder effects are less im- 
portant) than would be possible without it. 

Experimental evidence for the WC in 2DEG's has ac- 
cumulated over the last several years || . One probe that 
has produced much intriguing data is photoluminescence 
(PL), in which either a valence band hole Q] or a hole 
bound to an acceptor recombines with an electron in 
the 2DEG, producing a characteristic photon spectrum. 
A mean-field analysis [jj of the latter type of experiment 
showed that the PL spectrum has, in principle, charac- 
teristic signatures of the WC: a "Hofstadter butterfly" 
H spectrum for the case of weak interactions between 
the electrons and the hole, and a characteristic shift in 
the PL spectrum upon melting of the crystal. 

In this Letter, we go beyond the mean-field approxima- 
tion, to examine shakeup effects on the PL spectrum; i.e., 



we will examine how the collective mode spectrum of the 
WC (which, at long wavelengths, corresponds to the clas- 
sical phonon spectrum) , and the fact that some of these 
modes may be excited in the electron-hole recombination 
process, modify the results of the mean-field theory. We 
will consider in detail only the case of a localized hole 
Hf|. O ur method is purely quantum-mechanical, and 
treats both the tunneling electron and the other lattice 
electrons on the same footing. Furthermore, we employ 
a quantum treatment of the collective modes to realisti- 
cally account for contributions both from small and large 
wavevector collective excitations of the lattice. Since we 
are working in the strong magnetic field limit, we consider 
only excitations within the lowest Landau level (LLL). 
Our main results are: (1) Shakeup effects shift the main 
PL peak to higher energies than found in a mean-field 
treatment . (2) The Hofstadter spectrum is eliminated 
from the PL spectrum, even in the case of weak electron- 
hole interactions Q (although we will argue below that 
it survives in the itinerant hole case @). The sudden 
shift of the PL spectrum upon melting, by contrast, sur- 
vives even when shakeup is included. (3) Phonon side- 
bands appear that correspond to maxima in the phonon 
density of states (DOS); some (but not all) of these side- 
bands are results of van Hove singularities in the DOS, 
and so are characteristic of an ordered WC state. For the 
case of weak electron- hole interactions || , we do not see 
these sidebands in the liquid state, so that phonon satel- 
lites uniquely distinguish between a liquid and a solid 
state. Interestingly, for stronger electron-hole interac- 
tions, a shakeup satellite persists even above the melting 
temperature. We expect this sideband to lose oscillator 
strength relative to the main peak, either with increas- 
ing temperature or decreasing electron-hole interaction 
strength. The latter may be accomplished by examining 
PL from several samples with different acceptor - 2DEG 
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setback distances ||. 

Examples of our calculated PL spectra are shown in 
Fig. 1 for filling fractions v = 1/5 and v — 2/7, for elec- 
tron density N s = 6 x 10 10 cm -2 . Our hole is assumed 
to be strongly localized, and located 250A from the elec- 
tron plane. For the case of no electron-hole interaction 
[Fig. 1(a)], at low temperature, a well-defined shakeup 
peak may be seen approximately 2 meV below the main 
PL peak; a second very weak satellite is observed ap- 
proximately 3.5meV below the main peak. The origins 
of these peaks may be understood in terms of the phonon 
DOS, which is illustrated in Fig. 2. A van Hove singular- 
ity, arising from zone-edge phonons, appears as a strong 
double peak near 0.4meV. Two other peaks may be seen 
near 1.2 meV and 1.9 meV. There are weak sidebands 
associated with each of these peaks in the PL spectrum. 
The precise interpretation of these peaks is unclear; how- 
ever, it has been speculated that these represent vacancy- 
interstitial excitations |l0[ . We point out that it is crucial 
to use a fully quantum mechanical treatment of the col- 
lective excitations of the lattice to observe these higher 
order satellites; classical treatments of the phonons [p|JTl| 
do not produce these unusual excitations. 



to an upward shift of the main PL peak. Our calculation 
shows that the final states in which phonons are excited 
are not nearly so strongly renormalized by lattice relax- 
ation effects, leading to the increased splitting between 
the main PL peak and the sidebands. 

The inset to Fig. 1 illustrates the PL spectrum in 
the melted state both without and with an electron- 
hole interaction. In the former case, there is no phonon 
sideband present. This is necessarily so, because in the 
melted phase, the density is uniform, and there are no 
collective modes in the LLL ftofl . By contrast, when 
an electron-hole interaction is present, there is a non- 
uniform electron density near the hole, allowing some 
local collective modes to persist even above the melting 
temperature. With further increase in temperature, or 
increased setback between the hole and the 2DEG, the 
oscillator strength of this mode will significantly decrease. 

Fig. 1(c) illustrates our results for v = 2/7 in the 
absence of electron-hole interactions. As can be seen, 
except for a change in energy scale caused by changing 
the magnetic field, the lincshape is essentially identical 
to the case of v = 1/5. This contrasts sharply with the 
results found in the mean-field approximation M, where 



FIG. 1. PL spectra for (&)v = 1/5, T = 0, no elec- 
tron-hole interaction; (b)z/ = 1/5, T = 0, with electron-hole 
interaction; (c)u — 2/7, T = 0, no electron-hole interaction. 
Inset: PL spectra for (a) and (b), with T just above the melt- 
ing temperature. In order to distinguish different spectra, 
they were seperated by 250 units. 

It should also be noted that the splitting between the 
main PL peak and the first sideband is actually larger 
than the energy at which the van Hove singularity in 
the phonon DOS appears. The reason for this is that 
there is a strong self-energy renormalization due to the 
phonons in the main PL peak. Physically, this arises be- 
cause the final state of the crystal contains a vacancy, 
which is lowered in energy by a distortion of the lattice 
- i.e., by allowing the electrons surrounding the vacancy 
to relax inward. The self-energy shift accounts for this 
lowering in energy of the final state of the WC, and leads 



FIG. 2. Collective mode DOS for v = 1/5, T = 0. 

without electron- hole interactions, a filling v — p/q gen- 
erally yields p distinct lines for a localized hole. While 
the splittings are so small in that situation that they are 
difficult in practice to resolve, evidently shakeup effects 
wipe out this structure even in principle. We note that 
for the case of an itinerant hole, we expect these charac- 
teristic splittings to survive shakeup effects. The reason 
is that (neglecting excitonic effects), the PL spectrum is 
related to the product of Green's functions for the hole 
and the electron. For localized holes, only the latter has 
poles in the form of a Hofstadter spectrum at the mean- 
field level, which are wiped out by shakeup effects. How- 
ever, for itinerant holes, the hole Green's function also 
has poles of the Hofstadter form, which are unaffected 
by shakeup. Thus, one should in principle be able to 



2 



REVTgX 3.0 



Das Sarma Group Preprint, 1994 



identify this characteristic spectrum of bands and gaps that is unique to a WC in a magnetic field in itinerant hole 
experiments [l2| . 

We now present an outline of how these results were derived; a more detailed description will be presented in a 
future publication |l3[ ]. Previously 0, it was shown that if we consider a superlattice of holes, with n e electrons and 
one hole per unit cell, the PL intensity may quite generally be written in the form P(u>) oc lm[R(cu + i5)], where 
R(ui) = (71^/271-/0) ^ G i?(G,w)e _G z °/ 4 , and we have approximated the core-hole wavefunctions as delta-functions, 
n/j is the density of holes, the volume of the system, and the vectors G are the reciprocal lattice vectors of the 

/ \l/2 

superlattice. The magnetic length Iq — Ihc/eB) will be set to unity in the remainder of this paper. R(G,uj) is 

defined as R(G,u>) = (l/ff) Ylx e~ %GxX+tG * Gy l 2 Ru(X, X — G y ; u), where g is the Landau level degeneracy, X is the 
guiding center quantum number, and 

R lJ {X 1 ,X 2 ; icj n ) = I < T T a Xl {T)c t {r)c]a ] X2 > e l ^ T dr, 
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with a x creating an electron in state X and c\ creating a hole in the unit cell i, and /3 is the inverse temperature. 
Working in the lowest Landau level, the equation of motion for Rij(G,u>) may be written as 

J^(G,t)=< P (G,t) > 8 13 5{t) - e %(G,r) - n h J2 F(G')e iG ' xG / 2 - G ' 2 / 4 i^(G - G',r) 

G' 

~ E Kq)C^(-q, q + G)^ G ^/ 2 ^ 2 /2 _ _L £ 7(q)C^(-q, G)e-^^'\ (1) 

q#0 q 

In Eq.(|l|), w(q) and V(q) are the Fourier transforms of the electron-electron and electron-hole interactions, respectively, 
60 is the energy of the localized hole, and the sum over G' is only over reciprocal lattice vectors, while the sums over 
q are over all wavevectors. Ri specifies the position of the hole in the i th unit cell, and < p(G, r) > e~ G / 4 is 
the expectation value of a Fourier component of the electron density. While this quantity is independent of r in the 
groundstate, it will be convenient for later purposes to formally leave it as an argument of the density. The method 
for computing these Fourier components has been described elsewhere [[7 10 . The correlation function CV/(pi, P2) in 



the last two terms is just Fourier transformation of the following correlation function: 

C iS (X 1 X 2 ;X z X«r',T) = -9 < T T a x y)a X2 (T')a x ,(T)c l {T)c){Q)a Xi (0) > . (2) 

Eq.(Q) represents the first in an infinite series of equations relating an n particle Green's function to the n + 1 
particle Green's function |fl4|j . In the mean-field approximation, it is simplified by employing a Hartree-Fock (HF) 
decomposition of Eq.(||), which converts Eq.(jl]) into a self-consistent equation for i?^ 0. To include shakeup effects, 
we instead extend this hierarchy to one more level, writing down a self-consistent form for Cij which explicitly 
contains the collective mode excitations. To carry out this program, it is convenient to implicitly define a self-energy 
by substituting the last two terms in Eq.(§) by - J2 G > So dr'E(G, G'; r - r')%(G' ', r'). The HF approximation for 
the PL is equivalent to taking S(G, G'; r - r') by S HF (G, G')S(t - r'), where 

Z HF (G,G') = W(G - G') < p(G - G';r) > e tGxG '/ 2 + iVy(q) < p(-q,r) > s > (3) 

Z7T * — ' 



where W is the sum of the direct and exchange Coulomb potentials [|7|,|10[] . To generate a self-consistent equation for 
Cij, we take a functional derivative E3] of Eq.(|]) with respect to a spatially periodic external potential U(p, r'). In 
doing this, it must be noted that a term directly coupling the Green's function Rij to the external potential U must 
be added to Eq.(|l|), and the sums over reciprocal lattice vectors must be extended to all wavevectors, because U is 
not in general commensurate with the lattice. After an arduous calculation, the result expresses Cij in terms of the 
density-density correlation function x(Gi + q, G 2 + q; r) = XG!G2(q; T ) = — 9 < Tp{G\ + q, r)p{— G2 — q; 0) >. 
The Fourier transform of this response function contains poles at the collective mode frequencies of the system; it is 
thus through x that shakeup effects are introduced. We have computed this function numerically using a generalized 
random phase approximation p| which does not assume small displacements of the lattice electrons (as is necessary 
in classical approaches |^]), and thus gives a realistic dispersion relation for the collective modes across the entire 
Brillouin zone. Upon substituting the resulting Cy into Eq. (Q) , one obtains the self-energy S = Y, hf 5(t — t') + ST,, 
with 
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<5S(Gi, G2, T\ — T2) = — V / d 2 qe lqx ( Gl - G2 )/ 2 - lGl><G 'i/ 2+lG2><G ^ 2 iy(G' 1 +q)V7(G2+q) 

xF(Gi - Gi q, G 2 G' 2 - q; n - t 2 )xg;g^ (q; n - r 2 ) 

+ 7Al E / d 2 q e - (G ' 1+q)2/4 - (Gi+q)2/4 ^(G' 1+ q)F(G 2 +q) 

xF(Gi,G 2 ;n -r 2 )xG' 1 G' 2 (^n -r 2 ) (4) 

In Eq.(|4|), J BZ d 2 c[ represents an integral over wavevectors in the first Brillouin zone of the superlattice, and F is a 
generalized Green's function satisfying the equation of motion 

|-F(G, G lf r - n )= £ GGi <5(t - n ) - n fc J] ^(G V G ' xG / 2 - G ' 2 / 4 F(G - G', G x , r - n) 
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e F(G, Gi, t - ti) - 53 £ HF (G, G')F(G', G x , r - n), (5) 

G 



Eq.(||) may be solved using methods discussed previously RM. With this expression, we are now able to compute 
the PL intensity. We substitute T = T hf 8(t — t') + ST into Eq. ([!]), and Fourier transform this with respect to 
imaginary time. This means that a Fourier transform of Eq.(^) will be necessary, leading to frequency summations 
of the form (suppressing wavevector arguments) ^2 iuJ F(ioj n )x{<X, — iw„). To accomplish this, we represent \ as 
a sum over its collective mode poles E3j; the frequency sums may be then computed using standard methods p4[ . 
The computation of ST is clearly the bottleneck in this computation, since it requires two reciprocal lattice sums 
and an approximate sum over the continuous wavevector q. We have accomplished this using 469 q points in the 
first Brillouin zone, for one and three electrons per unit cell, which give very similar results. Finally, once we have 
computed ST, it is straightforward to substitute this into the frequency version of Eq. ([!]), obtain R(G,uj), and from 
there compute the PL spectrum. 

In summary, we have developed a method by which shakeup effects in the PL spectrum of a WC from localized 
holes may be computed, that treats the tunneling electron and the lattice electrons on an equal footing, and uses a 
fully quantum treatment of the collective modes that is realistic over the entire Brillouin zone. Our method is quite 
general, and should be applicable to other shakeup problems where quantum fluctuations are important. We find that 
the Hofstadter spectrum found in a mean-field analysis of this experiment is lost (although we expect it to survive in 
itinerant hole experiments), and is replaced by a series of sidebands due to creation of phonons and other collective 
excitations of the WC. These sidebands are a unique signature of the WC, and can in principle be used to distinguish 
between a liquid and crystal state of the electrons. We find that there is a sudden shift in the PL spectrum upon 
melting of the crystal. 
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